Vertical transport and domain formation in 
multiple quantum wells 

Andreas Wacker 

Mikroelektronik Centret, Danmarks Tekniske Universitet, DK-2800 Lyngby, 

Denmark 
email: wacker@mic.dtu.dk 



to appear in: 

Theory of Transport Properties of Semiconductor Nanostructures 

edited by Eckehard Scholl 
(Chapman and Hall, in preparation) 



In this book article effects related to the vertical transport in 
weakly coupled multiple quantum wells are reviewed. A self-contained 
microscopical model for the calculation of the well-to-well currents 
without any fittings parameters is presented. The model exhibits the 
well-known peaks in the current-field relation in quantitative agree- 
ment with experiments. This local current-field relation is used as an 
input for the calculation of the transport in the extended structure 
consisting of many periods. Here both the formation of stationary 
field domains as well as self-sustained current oscillations are found 
in good agreement with experimental data. The underlying physics 
of these nonlinear phenomena is discussed in detail. 
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Introduction 



Today's growth techniques allow the construction of semiconductor structures 
where layers of different semiconductors (exhibiting similar lattice constants) 
can be grown on each other with the interface being well defined within one 
atomic monolayer. If two such layers alternate periodically one obtains a pe- 
riodic structure with an artificial period d in the growth direction (which is 
defined to be the ^-direction) . This leads to spatial variations in the conduction 
and valence band of the material as sketched in Fig. |j} Considering only con- 
duction band states in the following, the region with lower conduction band is 
called 'well' and the region with higher conduction band 'barrier'. An extended 
discussion of various aspects of such structures can be found in Ref. Q] . In this 
chapter the vertical electronic transport (i.e., in z-direction) is considered for 
such structures. 

If we consider the electronic properties of such a structure, there are two dif- 
ferent approaches. In the first approach one may consider the full structure as 
an artificial lattice. The energy spectrum can be calculated analogously to the 
Kronig-Penney model (which is discussed in almost all solid-state physics text- 
books) resulting in the appearance of energy bands and energy gaps as sketched 
in Fig. |l} Due to this analogy with the atomic lattice of lattice constant 
such semiconductor structures are often called superlattices. The correspond- 
ing eigenfunctions are the usual Bloch functions which extend over the whole 
superlattice. The Bloch functions are labeled with the Bloch vector q which is 
restricted to the Brillouin zone — Tr/d < q < ir/d. This range is much smaller 
than the Brillouin zone —n/a^ < q < t^/cll of the atomic lattice as d > a^. 
Therefore these new bands are called minibands. If an electric field is applied, 
the Bloch functions are accelerated and the transport can be treated analogously 
to the usual transport in bulk systems Q . Due to the short Brillouin zone many 
special effects like Bloch oscillations can be found here which are not accessible 
in bulk systems as discussed in the preceding chapter. 

In the second approach one may consider the wells to be isolated from each 
other as a first approximation and calculate the eigenstates within each well. 
This yields a sequence of energy levels and localized wave functions for each 
well. Such a structure is usually called a multiple quantum well. Of course 
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Figure 1: Sketch of the conduction band E c (z) with minibands v = a,b 
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this approach only makes sense if the coupling between different wells is weak 
and can be calculated perturbatively. Thinking in terms of Fermi's golden rule, 
there will be transitions between the levels in different wells yielding sequential 
tunnelling. Due to energy conservation this tunnelling takes only place if the 
energy levels align and strong resonances || are likely to occur if the relative 
height of the levels in different wells is varied by a voltage applied in the z- 
direction. 

The width of these resonances is related to scattering processes of the elec- 
trons inside the wells. In Ref. Q this feature was used to determine the electron- 
electron scattering rate from the transport between two weakly coupled quantum 
wells. This shows that a modelling of scattering processes is necessary for a cal- 
culation of the current in such structures and opens the possibility for checking 
theoretical concepts as well as investigating the significance of various scattering 
mechanisms by comparison with experimental data. 

These resonances yield strong nonlincarities in the local current-field rela- 
tion. If we now have a long periodic sequence of these coupled quantum wells the 
full system is a good example of an extended nonlinear system. This leads to the 
appearance of complicated current-voltage characteristics exhibiting an almost 
periodic sequence of branches due to the formation of electric field domains as 
well as to self-sustained oscillations due to moving domains. As the periodicity 
of the structure is an important feature for the translational invariance, these 
structures are frequently called superlattices for weak coupling as well. If the 
translational invariance is broken due to an insufficient control of the growth 
conditions the branches lose their periodic structure, and information about 
the actual structure can be obtained from measurements of the current- voltage 
characteristics. 

The first section of this chapter is devoted to the classification of the different 
regimes depending on the coupling between the wells. In the second section a 
transport model for weakly coupled quantum wells is presented which allows a 
microscopic calculation of the current without any fitting parameters. In the 
third chapter this model is extended in order to describe domain formation. The 
fourth section is on the influence of deviations from the periodic structure of 
the superlattice and the fifth section refers to self-sustained oscillations. Finally, 
details of the calculations will be given in the sixth section. 

1 The different transport regimes 

In this section the precise meaning of weakly and strongly coupled quantum 
wells is discussed. 

Assuming ideal interfaces the semiconductor structure is translational in- 
variant within the x and y direction perpendicular to the growth direction. 
Therefore the x,y dependence can be separated by the ansatz e-- where k 
and r are vectors within the two-dimensional (x, y) plane. (The effects due to 
the atomic periodicity of the lattice are treated within the envelope function 
formalism H assuming a parabolic subband.) Within the ^-direction we have 
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Figure 2: Conduction band E c (z) together with the Wannier functions calcu- 
lated for the two lowest subbands. 



an artificial period d leading to energy states E" characterized by the miniband 
index v and the quasi- wave vectors — ir/d < q < ir/d. The energies as well 
as the wave- functions <fiq{z) can be calculated numerically analogously to the 
well-known Kronig-Penney model. For a given miniband v the energy E v q has 

a mean value E v = d/{2n) f^^dqE" and varies within the miniband width 
A" = Max g (££) - Mirig(^) as sketched in Fig. [| 

From the Bloch functions (fiq(z) the Wannier functions ty"(z — nd) defined 

by 



^ v {z - nd) = 




can be constructed. They are real and localized in well n for an appropriate 
choice of the complex phase for ip v q (z) Q. An example is shown in Fig. ||. Using 
the Fourier expansion 

oo 

E»=E» + J22T»cos(hdq) (2) 

h=l 

one obtains the following Hamiltonian in second quantization 

oo 

H = £ + ]T (T£<&„< + T^J_ h cl) (3) 

n,v h=l 

for the creation d£ and annihilation c v n operators of the Wannier functions 
^ u (z — nd). As the Wannier functions are linear combinations of Bloch functions 
with different energies, they are no stationary states. Neglecting terms with 
h > 1 the time evolution of the annihilation operators is given by 

;/4< = rr (<+!+<_!) . (4) 
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For the initial condition c^(t = 0) = S n ,oCQ this has the solution 

<(t)=rV n (—!-*) eg (5) 

where J n (x) is the Bessel function of order n. Since YlnL-oo n2 {Jn(x)) 2 = x 2 /2 
the average extension n a (t) as a function of time is found to be 

V2T X 



n a (t) = ^) = ^_±t. (6) 

This is obviously a coherent process. In a real semiconductor structure there 
will be scattering (characterized by a scattering time r sc ) destroying the co- 
herent evolution^ given by Eq. (Q). If n a (r sc ) -C 1 holds, the phase coherence 
is completely lost during the spread and the states of adjacent wells will not 
maintain a fixed phase relation. In this case the transitions can be described 
by sequential tunnelling between the Wannier states with loss of phase between 
different tunnelling events. Defining T = h/r sc one therefore may use the picture 
of 



sequential tunnelling for T\ -C T/y/2 . (7) 



If, on the other hand, the extended Bloch functions shall be a reasonable 
basis set, the phase coherence has to be maintained over a larger number of 
quantum wells. This gives n a (r sc ) > 1 as a first necessary condition. But a 
further complication arises if an electric field F is applied in the z direction. As 
discussed in the preceding chapter of this book the eigenvalues of the Hamilto- 
nian take discrete values E v — neFd which is called the Wannier-Stark (WS) 
ladder. The corresponding eigenfunctions 4>^ s (z — nd) are localized around well 
n. Within the one band limit they are given by 

-LU»( z -nd). (8) 

(See, e.g., Ref. jjj, where also the coupling between different bands is discussed.) 
Then the spatial extension of the WS-states is given by 



n ws = ^) = ^. (9) 

If now nws becomes small, the electrons are localized and the use of extended 
Bloch functions does not make any sense. This provides the second condition 
^ws ^ 1 for the Bloch functions to be an appropriate basis set of wave functions. 
Thus it only makes sense to speak about 



minibands for T x > T/V2 and T x » eFd^/2 . (10) 



1 This scattering can be a phase breaking process such as scattering with phonons. But also 
elastic scattering processes at impurities or interface roughness will contribute here because 
they destroy the translational invariance of the superlattice, which is both the justification 
for the existence of the Bloch functions and a WS state. 
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Figure 3: Sketch of the regimes where different approaches are valid 



As a third basis set the Wannier-Stark functions <Pws( z ~ nc ma y ^ e use d 
if an electric field is applied. They make sense if coherence is maintained within 
their extension, i.e., if ?i\vs < n a (T sc ). Therefore the transport can be described 

by 



WS hopping for eFd > V . (11) 



These different regimes of validity are sketched in Fig. || where the condition 
a <C b has been translated to a < b/2 for illustrative purpose. 

How can the electric transport be described if an electric field F is applied 
in the z direction? For F = the system should be in equilibrium and no 
current / flows. For small fields the current should increase with F. (This 
may be calculated by linear response which is given by the Kubo formula for 
a quantum system. See, e.g., chapter 5, where the application of the Kubo 
formula has been described for the Coulomb drag.) Now it is an important 
feature of both superlattices and multiple quantum wells that for larger fields 
there is a range with negative differential conductivity dl/dF < (NDC). This 
can be easily understood within the Wannier-Stark picture. As the WS-states 
are orthogonal, transitions between them have to be caused by scattering events. 
Now the matrix element for a given scattering element depends on the spatial 
overlap of the eigenfunctions. With increasing field the Wannier-Stark states 
become more and more localized as can be seen from Eq. (^). Therefore the 
matrix element for the scattering is decreased and the hopping rate is diminished 
H ||. As the WS approach is not justified for low electric fields (see Fig. 0) an 
Ohmic behaviour for low fields could not be observed within this model M. 

In the validity range of ([!(]) the basis of Bloch functions may be used. Due to 
an electric field the Bloch waves are accelerated according to Hq — eF. Like in 
the standard transport theory for bulk systems a positive conductivity is found 
for low electric fields. But for higher fields the Bloch vector q may run through 
the whole Brillouin zone 2ir/d if eFT sc /H > 2ir/d holds. Then the electrons 
perform a periodic motion in both q and z space which is called the Bloch 
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oscillation. A quantitative analysis within various approaches j^, 10, [IT], p"2| 
reveals 

eFd 

(eFd)*+r*- (12) 
This yields positive differential conductivity for eFd < T and negative differen- 
tial conductivity for larger fields which is confirmed by experimental data (see, 
e.g., [|l3| Q). For large F there is / ~ l/F. This behaviour is also found from 
the WS-hopping model for eFd < A = 4Ti § which can be easily understood 
from Fig. |§| as this condition includes the region where both approaches are 
valid. This is an explicit example of the more general equivalence between the 
Wannier-Stark picture and the Bloch oscillation picture discussed in chapter 9 
of this book. 

In case of weakly coupled multiple quantum wells T\ <C T/\/2 NDC has 
been observed experimentally as well jl5fl . Here we can consider the electronic 
states to be the levels of the single quantum wells which will exhibit a certain 
width r due to scattering. The tunnelling between the wells is caused by the 
coupling T\. Due to energy conservation a significant current between adjacent 
wells may flow if the levels are aligned within accuracy T. Therefore NDC is 
likely to appear if the disalignment eFd becomes larger than T |l6|, [l7| . This 
will be modelled quantitatively in the next section. 

For the white region of Fig. || between these two limits the situation is more 
complicated. Here Laikhtman and Miller |l8| were able to show that Eq. ( [i"2| ) 
holds provided the electron temperature is much larger than T\ and eFd. 

In an extended system the existence of an NDC region typically causes os- 
cillatory behaviour due to travelling field domains (see, e.g., For the mini- 
band regime this has been studied theoretically quite early in Refs. |^0[ ^lj 
and recently also experimentally |23[ p3. Note that these travelling domain 



oscillations are self-sustained in Ref. 124] while for the usual Bloch oscillations 
only transient oscillations are found (see, e.g., Ref. (^5|). The same type of 
oscillations has been recently found in the regime of sequential tunnelling both 
experimentally |2(| |2^] and by numerical simulations |p9| , |3Cf| . An analytic 
treatment of the oscillation mode is given in Refs. ^]. Under certain con- 
ditions chaotic oscillations occur as well ||| [34) . 

For weakly coupled multiple quantum wells an additional scenario occurs, 
which has been extensively studied in the last decade: As the electronic states 
are localized within single wells the domain boundary (which are charge ac- 
cumulation layers) can be trapped in a single quantum well and a stationary 
stable domain is formed. As the domain boundary may be located in any well 
a periodic sequence of branches appears in the current-voltage characteristic. 
This effect has been observed experimentally by many groups ]35|]- [^9|. The- 
oretically such effects can be studied by combining rate equations between the 
wells and Poisson's equation. Such an approach has already been performed 
in Refs. (5(], |5l| . Nevertheless, to my knowledge, the full current- voltage char- 
acteristic exhibiting the domain structure could first be resolved qualitatively 
in Ref. [p2| for slim superlattices exhibiting a few electrons per well and in 



Refs. p3l for two-dimensional wells. A full quantitative calculation has been 
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presented recently |35| [5(| for different multiple quantum wells yielding quanti- 
tative agreement with experimental data without using any fitting parameters. 
This model will be presented in the next sections. 



2 Modelling of the transport between weakly 
coupled quantum wells 

In this section I want to show how the currents between the wells can be cal- 
culated from a microscopic model. Furthermore I present some simplifications 
in order to obtain estimates for the current. The numerical calculations are 
performed with the data of the sample used in Refs. [^6[ [to| exhibiting TV = 40 
GaAs wells of width w — 9 nm between 41 AlAs barriers of width 6 = 4 nm. 
The wells are n-doped with a doping density of Nj) — 1.5 x 10 11 /cm 2 per well. 
Details of the specific calculations outlined in this section are given in section 

As discussed in the first section, for weakly coupled quantum wells the prod- 
ucts of Wannier functions ^l"{z — nd), localized in well n, and plane waves e — - 
form a reasonable basis set of wave functions. Now we restrict ourselves to the 
lowest two levels denoted by v — a, b. The respective annihilation operators are 
denoted by a n = c° and b n — c b n . Furthermore we add the contribution —eFz 
due to a homogeneous electric field F in the Hamiltonian (||) where e < is the 
charge of the electron. Restricting ourselves to coupling between neighbouring 
wells (h = 1) we obtain the Hamiltonian H = H + H\ + H 2 : 



H o BS = ^[(£ Q + £ fe -eFdn)4(fcW(fc) 



n.k 



+ (E b + E k - eFdn)bl(k)b n (k)] (13) 
H{ cs = Y,[ T >n + Ma n {k) + Ttbi +1 {k)b n {k) 

-eFRfal +1 (k)b n (k) - eFR\ a bl +l {k)a n {k) + /i.e.] (14) 



with the parabolic dispersion Eu = h 2 k 2 /(2m w ) {m w is the effective mass in the 
well) and the couplings R% v = J dz^ v (z — hd)z^/^(z). All energies E are given 
with respect to the bottom of the quantum well. The values of the coefficients 
are presented in Table [l]. The term H 2 can be incorporated into the one electron 
states by diagonalizing H + H 2 . This leads to renormalized field-dependent 
coefficients in H and Hi (which are used in the following) but does not change 
the structure of the problem for a homogeneous electric field. 

H rcB is only considering ^-conserving processes reflecting an ideal structure 
with translational invariance in the r-plane and neglecting any many-particle 
processes. Nevertheless, there are non fc-conserving processes as well which may 
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E a = 47.1 mcV 


Tf = -0.0201 mcV 


R b a = -0.14M 


E b = 176.6 meV 


= 0.0776 mcV 


R\ a = 2.66 x 10" 4 d 



Table 1: Calculated level energies and transition elements for Eqs. (|13||l 



result from scattering at impurities or interface roughness, e.g.. These give two 
further contributions to the Hamiltonian: 

^scatter 1 £ [[/-(^t (fc + p)a n (k) + U$> {^)b\{k + p)b n (k) 



h_.jl 



(16) 



+U*?(p)bl(k +p)a n (k) + US b (p)ai(k+p)b n (k)] 
is the contribution due to scattering within the well and 

£ fa tter = 1 [ur(p)ai +1 (k+p)a n (k) + Ut b {p)bi +1 {k+p)b n {k) 

k,p 

+U b 1 a (p)bl +1 (k+p)a n (k_) + U1 h {p)ai +l (k+p)b n {k) + h.c] (17) 

refers to interwell scattering, where the restriction to neighboured wells is made. 

Within the assumptions of local thermal equilibrium in each well and weak 
coupling between the wells the current from level v in well n to level fj, in well 
n + 1 is given by the following expression^ (see section 9.3 of Ref. p8|): 



(18) 



< +1 (k', E + eFd) [n F (E - E%) - n F (E - < +1 + eFd)} 



Here E^ is the chemical potential in well n which is measured with respect to 
the bottom of the quantum well. n_p(x) = (\ + e x ^ kBT '})~ 1 is the Fermi function 
and T e is the electron temperature. A^(k,E) denotes the spectral function for 
the state k of the subband v in well number n. It represents the weight of the 
free particle state k contributing to the energy E. Then the total density of 
states Pn(E) in subband v is given by 



Pn{E) 



2ttA ^ 

k 



Kik.E) 



(19) 



where the factor 2 reflects the spin degeneracy and A denotes the sample area. 
If no scattering is present, the state k has a fixed energy E v + Ek, and the 



2 An important point of the derivation is the assumption of uncorrelated scattering in 
different wells. For example, this is the fact if the electrons are dominantly scattered by the 
impurities localized in the same well. If the scattering occurs at identical impurities for the 
electrons in well a and in well n + 1 correlation effects occur which may essentially change the 

result, sec Refs. iplna. 
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E k [meV] E [meV] 

Figure 4: (a) Spectral function A a (Ek,E) of the lowest level for two differ- 
ent energies E calculated for impurity scattering and interface roughness (see 
section ||). (b) Density of states in units of the free-particle density of states 
(dashed line) where we have also indicated the value of the Fermi-energy at zero 
temperature. 



spectral function becomes a 8- function A^(k, E) = 2n8(E — E v — Ek). In the 
continuum limit (J2k ~ * A/(2ir) 2 f d 2 k) Eq. ( pj| ) gives then the 2-dimensional 
density of states p v n {E) = p Q 9{E - E v ) with p = 

If scattering is present due to H g catter the states k are no longer eigenstates 
of the total Hamiltonian. This can be taken into account by calculating the self- 
energy for the given scattering within standard theory (see, e.g., (58)). Assuming 
equilibrium the spectral function is then related to the retarded self-energy 
T,» let (k,E) via 

-2ImE" et (fc,£;) , . 

A v n (k,E) = " v ~' ; j ■ 20 

{E-E v -E k - Re££ ret (fc, E)f + (Im££ rct (/c, E)) 

Frequently, the self-energy is taken to be constant for simplicity setting S^ rct (k, E) 
W% — Then the spectral function is a Lorentzian with a full width at 

half maximum T. 

Here the self-energies are calculated from basic scattering processes at im- 
purities and interface roughness without using any fitting parameters. This 
calculation is presented in section p|. The calculated spectral function for differ- 
ent energies E are shown in Fig. RFa) for illustration. One can clearly see that 
they exhibit a maximum close to Ek — E. From the width we may estimate 
that r w 10 meV holds for E = 2 meV and T w 6 meV for E = 10 meV. 
Fig. ||(b) shows the total density of states. 

While the full derivation is slightly tedious [Q formula (|l^) can be motivated 
quite easily: In the long-time limit energy has to be conserved during transitions 
caused by the time-independent interwell couplings H^ +1 -^ k , nk . Therefore we 
have to consider tunnelling processes for a certain energy E and integrate over 
E afterwards. The factor [np(E — E^) — np(E + eFd— E^ +1 )] takes into account 
the thermal occupation at the given energy in both wells. The free particle state 
k has a weight A v n [k, E)/(2n) in well n. Its transition probability to the state 
fc' in well n + 1 is given by 2Tr\H^ n '^_ 1 ^ k , nk \ 2 ' Jh (Fermi's golden rule). The final 
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state has a weight A^ +1 (lc , E + eFd)/{2ir) at the given energy. Obviously one 
has to sum over all free particle states k, k. Finally, the factor 2 is due to the 
spin degeneracy. 

2.1 Resonant Transitions 

Let us first investigate the current due to the transition elements from H[ es in 
Eq. (^). T hey conserve the momentum k and therefore the kinetic energy Ek- 
Using Eq. (^9|) one can rewrite Eq. (|lj) as 

I jrtesn,v,2 roc 

CT£T ^ A \ J JEp» n {E)(A» n+l ){E,eFd) 

[n F (E-E^)~n F (E~E^ +1 +eFd)] (21) 

with the average 

J2 k K(k, E)A»{k, E + eFd) 
(An +1 )(E,eFd) = "'- E ^;-) (22) 

If the spectral functions are 5-functions (A^ +1 )(E, eFd) = 2ir5(E l/ + eFd — E^) 
holds. Therefore tunnelling only takes place if the levels are exactly aligned. In 
order to estimate the effect of broadening one may assume constant self-energies 
E^ rct = —iT v /2. Performing the continuum limit and assuming E — E v ^> 
as well as E + eFd - E^ > one finds 

T v + V 1 

(K+i ) {E, eFd) = {eFd + Ev _ E , )2 + {TV + T , ) 2 /A ■ ( 23 ) 

This expression has a peak at the resonance E 11 = eFd + E v and a full width at 
half maximum of {T u + r M ). Even if the conditions stated above are not fulfilled, 
the result is typically similar. 

For tunnelling between the lowest levels a —> a, equal densities E^ — 
E^ +1 , a constant density of states p a and not too high temperatures and fields 
k B T, eFd < El - E a one finds 

r - - cA A T ? )2 2TaeFd (24) 
l n ^ n+1 -eAp h {eFd)2 + {Ta)2 ■ (24) 

Therefore the current has a maximum at eF m&x d = T a where it takes the value 
Im&x — eAp a (Tf) 2 /h. Note that this value neither depends on T nor on the 
Fermi level. Estimating T « 8 meV from the spectral functions shown in Fig. Q 
and taking p a = po yields 

eF max d sa 8meV / max « 0.27mA . (25) 

for the a — > a transition using Tf = —0.0201 meV from Table [l[ 

With the translations 2Tf -> E x , E a -> E , and h/T a -> r Eq. (|J) is 
identical to equation (12) f\ of Ref. |l0) for the case E F > E + E x . There 

3 Note that there is a factor r missing in the numerator due to a misprint. 
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semiclassical transport in a miniband for the strong coupling limit Tf 3> T a 
was regarded while here the limit of weakly coupled quantum wells Tf <C T a 
is considered. Note that in both derivations the case (E F — E a ) > T a ,2T 1 a is 
considered. 

The transitions from the lowest level to the excited level a — ^ b will become 
important if eFd w E b - E a . Assuming that k B T,E F +1 < E b - E a we find 
np(E — E F +1 + eFd) w in this field range. Using the approximation ( p3| ) we 
find close to the resonance 

I TrTesa,b\2 pa , p6 

ra^b .a ]££]__! 1 +i ( 2 a\ 

h ( eFd + E a _ £b)2 + (p, +I*)2/ 4 " ^ 

Note that this current is proportional to the density of carriers = J dEnpiE— 
E^)p a (E) in well n while this was not the case for the a — > a peak in Eq. (e3). 



2.2 Nonresonant Current 

The transition elements from //scatter ^ Q nQ ^ conserve momentum and therefore 
the kinetic energy Ek is changed during the transition. Thus these transitions 
are not as sensitive to the alignment of the levels as the resonant transition 
discussed above. They yield a background current which may dominate the 
current between the resonances. Here only nonresonant transitions via interface 
roughness are considered. The explicit expressions used are given in section |[ 
As the actual shape of the interface is not known and may strongly vary for 
different wafers it has been parameterized by a reasonable set of parameters in 
order to show the magnitude of the effect. Therefore the reader has to keep in 
mind that the nonresonant currents can easily vary by a factor of three or even 
more in the following calculations. 



2.3 Calculation of the current 



In Eq. (18) the current depends on the electric field and on the Fermi energies 
in both wells. The Fermi energies Ef can easily be calculated from the two- 
dimensional density of carriers rii in the well using the relation 

m= f dEn F (E - Ef) pP{E) ■ (27) 
J n 

Finally, the total current from well i to well i + 1 is the sum of the contributions 
between the different levels. 

h^i+i = £ /n+i = HP, ni,n i+1 ) (28) 

The result is presented in Fig. |B| for n, = = Nd- Neglecting the nonreso- 
nant transitions (dashed line) there is a first maximum at eFd = 9 meV with 
a current of 0.277 mA which is in good agreement with the estimation pa) for 
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Figure 5: Current Ii—>i+\ for = rij+i = Njj. The dashed line gives the 
current from resonant transition, while the full line gives then sum of resonant 
and nonresonant transitions due to interface roughness. 



the a — > a transitions. A second maximum occurs at eFd = 140 meV with 
/ = 1.75 mA. In comparison to this the experimental data (see Fig. 6 of |Q) 
exhibit a first maximum of / w 0.076 mA and a second maximum of / ~ 1.45 
mA. While there is good quantitative agreement for the height of second maxi- 
mum, the calculated first maximum seems to be too large by a factor of 4. This 
inconsistency will be resolved in the next chapter where the formation of field 
domains is considered. In contrast to the currents, the position of the maxima 
is much more difficult to compare as a part of the voltage may drop outside 
the superlattice. The nonresonant currents are not very large compared to the 
currents at the maxima but can dominate the total current between the maxima 
as can be seen from the full line in Fig. |^. 

Note that the experimental data of Ref. [Q indicate that the first excited 
level is unoccupied close to the onset of the a — > b resonance indicating that the 
current of approximately 0.1 mA is not carried by the a — > b transitions there. 
This may be attributed to stronger nonresonant transition in this field range. 

An important feature of the model presented here is the fact that only the 
nominal sample parameters are involved in the calculation of the currents. The 
resulting currents are in good quantitative agreement with the experimental 
data of Refs. Q. The same model has also been applied to the samples 
used in Ref. |39| where quantitative agreement could be obtained assuming a 
smaller barrier width |5f| . While these two highly-doped samples exhibit a den- 
sity of states which resembles the free-electron density of states with a smoothed 
onset (see Fig. |J), the density of states is much more complicated for low-doped 
samples due to the presence of impurity bands |^] . Within the single-site ap- 
proximation for impurity scattering used here (see section |^) these effects are 
included in the spectral function. The calculation yields a strong dependence 
of the current on temperature in good agreement with experimental data |60| . 
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Figure 6: Sketch of the fields and densities for domain formation. 



Excellent quantitative agreement has also been found when the sample was irra- 
diated by a strong terahertz field from a free-electron laser source . Therefore 
we may conclude that the formalism described here allows the quantitative cal- 
culation of the currents in weakly-coupled multiple quantum wells for a wide 
range of samples without applying any fitting parameters. 

3 Formation of field domains 

Now we want to consider the full structure consisting of N wells numbered 
i = 1, ...N and N + 1 barriers. Assuming that the electric field is constant 
over the structure the total voltage is given by U ~ (N + l)Fd. Fig. || shows 
that Ii-ti+i(F, Nd, Nd) exhibits a region of negative differential conductivity 
(NDC) between eF max d — 9.7 meV and eF min d = 90 meV. If the electric field is 
within this NDC-region an instability is likely to occur because a spontaneously 
formed charge accumulation will increase in time instead of decreasing. This 
is a common situation for samples exhibiting an N-shaped local current-field 
relation like the Gunn diode (see, e.g., jl^] and references cited therein). In 
order to include such effects charge accumulation has to be allowed for. This 
means that the carrier density rii in well i may deviate from the doping No per 
period. Integrating Poisson's equation over one period yields 

e r eo(Fi - Ft-i) = e(m - N D ) fori = l,...N (29) 

where Fi is the electric field in the middle of the barrier between wells i and 
i + Here a constant dielectric permeability e r — 13 is assumed for simplicity. 

The notation as well as a typical potential profile is sketched in Fig. ^| for 
a superlattice with N = 7. Within the first wells there is 1%% — rii — n% = ./Vp 
so that eFo = eF\ = eF^ = eF%. The electric field is low so that the current 
is dominated by the a — * a transitions. At well 4 there is charge accumulation 
71.4 > No. Therefore we have eF± > eFs and the electric field is large after well 

4 so that the current is dominated by the a — > b resonance. Here we assume 
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that the electrons relax fast from the upper levels into the lower level which is 
the only level to be populated in thermal equilibrium for level separations larger 
than the Fermi energy.^] 

In order to calculate the current between the wells the question arises which 
electric field should be used if an inhomogeneous situation is considered. As the 
voltage drop between the wells i and i + 1 can be approximated by eFid it is 
natural to use the electric field Fi in the argument of Eq. ( p8| ) . Nevertheless it 
has to be stated that this is an approximation done for simplicity and deviations 
may occur for inhomogeneous field profiles. 

Now the temporal evolution of the densities within the wells is given by the 
continuity equation 

eA— = I^_x)-^% — h-*{%+\) ^ I(F l - 1 ,n i - 1 ,n l ) - I(Fi,ni,n i+1 ) (30) 

for i = 1, . . . N. In order to obtain a complete set of equations for the temporal 
evolution of the densities and fields we have to add two more features. At first 
the voltage condition is now given by 

N 

U = dzF(z) m22dFi + U c (31) 

where U c represents the voltage drop outside the superlattice, which is neglected 
in the following in order to concentrate on the features of the pure superlattice. 
At second the currents Jo-a an d ijv->(jv+i) across the first and the last barrier 
of the structure, respectively, have to be specified. For simplicity one may use 
the expression ( pcf ) with appropriate effective densities 

no = n (F , N D ,ni) and n N+ i = n N+1 (F N ,N D ,n N ) . (32) 

In mathematical terms the functions no(Fo, Ni),ni) and un + i(Fn, Nd, n^) 
then represent the boundary conditions of the model. Eqs. (p^-p2[) form a 
complete set of equations for calculating the densities, fields, and currents as 
a function of time for fixed bias voltage U and given initial conditions niitg) 
(i=l, ; .N). _ _ 

Using a different approach for evaluating the current function I{Fi, rii, rii + i) 
(which is essentially based on the broadening due to the tunnelling time for 
the a — > b resonance and miniband conduction for the a — ► a resonance) and 
the boundary conditions uq(Fq, Njj,ni) — n%, iin+i(Fn, Nd, njsr) such a 

model has been used in Ref. |5j| to calculate the current-voltage characteristics 
under domain formation. The model of Ref. |54| uses the boundary condition 
rio(i 7 b, Nd, ni) = Njj and the current expression 

= eAn l v(F l ) (33) 



furthermore the intersubband relaxation must be fast enough to guarantee thermal equi- 
librium between the levels. If the level separation is larger tha n th e optical phonon energy 
36 meV, the intersubband relaxation time is of the order of 1 ps|6l[ , which is typically much 
faster than the tunnelling times. But for wide quantum wells the level separation becomes 
small and nonequilibrium effects are found|62|. 
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Figure 7: Calculated current-voltage characteristic exhibiting field domains for 
voltage sweep-up for the boundary conditions no = 3Njj and n^+i = 3Njj with 
(a) and without (b) nonresonant transitions due to interface roughness. The 
inset gives the field distribution at U = 2.5 V where the domain boundary is 
located at the 19th well. 

where v(F) is a phenomenological tunnelling rate yielding an N-shaped charac- 
teristic as shown in Fig. |^. This simplification allows the analytical construc- 
tion of the full current-voltage characteristic |54|, |(|| . Eq. ([33]) is motivated by 
Eq. ( |26| ) for the a — > b transitions and higher resonances as well. Nevertheless, is 
seems to be questionable close to the a — > a maximum as the maximum current 
of Eq. ( pi| ) is independent of rij for equal carrier densities in adjacent wells. The 
reason for this deviation is the existence of backward currents for low fields which 
are taken into account by the factor \jif{E — E%) — uf{E — E^ +1 + eFd)] in 
Eq. @. If on the other hand eFd^> E? +1 - E a + k B T e so that n F (E - E% +1 + 
eFd) w for all relevant energies E, Eq. ( |3^ ) is a reasonable approximation. 

3.1 Numerical results 

In a first simulation the boundary conditions uq (F, Njj ,n{) = n n+i (F, No , tin) = 
3Nd are used. The stationary stable states are determined by simulating 
Eqs. (|2S|-^2|) until a stationary state is reached for given U. Increasing or 
decreasing U afterwards simulates a sweep-up or sweep-down of the voltage, 
respectively. The result is shown in Fig. 0(a). We find that the characteristic 



16 



consists of 40 branches, equal to the number of quantum wells. The maximum 
current within the branches is about 0.09 mA in good agreement with the ex- 
perimental situation [l2| where around 0.06 mA is observed. This current is 
significantly lower than the first maximum o{I(F,N D ,N D ) (see Fig. |) resolv- 
ing the discrepancy discussed at the end of the last section. Note that another 
superlattice with the same specification p4| exhibits maximum currents of the 
branches of 0.14 mA. This indicates that the current is quite sensitive to vari- 
ations in the sample. The calculation exhibits similar effects. If we ignore the 
nonrcsonant transitions the maxima of the branches are found at 0.06 mA as 
shown in Fig. 0(b), although the height of the maxima in the homogeneous char- 
acteristics (see Fig. 0) is almost identical. Thus the extension of the branches 
are very sensitive to the quality of the interfaces (which causes nonresonant 
transition via interface roughness). Another uncertainty may be the actual bar- 
rier thickness as a variation of one monolayer changes the matrix elements H\ 
by a factor of 1.4 and therefore the current (which is proportional Hf i+l ) by a 
factor of 2. 

The slope of the branches varies between 1.1 mS and 0.55 mS for low and high 
voltage in Fig. which is significantly larger than the experimental slopes (see 
[fi"2"|) varying between 0.133 mS and 0.064 mS. This lower slope is responsible for 
the stronger overlap between the branches yielding pronounced multistability 
[ Ell . The discrepancy may be resolved by assuming an appropriate contact 
voltage U c in Eq. (^lj) which will depend on I and the fields at the boundaries. 
Another possibility might be that further nonresonant transitions alter the shape 
of the homogeneous characteristics of Fig. 0. 

The field distribution shown in the inset of Fig. 0(a) is in good agreement 
with cathodoluminescence measurements |35| stating that the field distribution 
consists of one low field and one high field domain where the high field domain 
is located at the anode. These measurements where performed at the same 
sample but in the range between the a — > b and the a — > c resonance. 

Using different boundary conditions almost identical domain branches are 
observed as shown in Fig. 0. The main difference occurs at the first and 
last branch which may be changed significantly. Comparing Fig. 0(a,b,c) and 
Fig. ^|(a) shows that the first branch is dominated by the boundary condition 
hn + i(Fq, Nd, ni) and the last branch by no(F]\[, Np, n^)- 

3.2 General aspects of domain formation 

Now a qualitative explanation of the domain formation shall be given. First 
regard Fig. H. For currents h from the interval [I m m, -fmax] where 7 ma x = 0.30 
mA is the maximum current for low fields and / m i n = 7.8/zA is the minimum 
current for medium fields there are 3 intersections Fj < Fjj < Fjji of the I(eFd) 
curve with Iq. Unlike Fjj the fields Fj and Fiji are in the range of positive 
differential conductivity (PDC). Therefore charge fluctuations are damped out 
in spatial regions where the electric field takes the values Fj or Fjjj. The field 
distribution for U = 2.5 V in the inset of Fig. 0(a) now shows that the electric 
field takes the value Fi ss Fi for i < 17 and Fi ss Fjjj for i > 20. The small 
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Figure 8: Calculated current-voltage characteristic for voltage sweep- up for the 
boundary conditions no = I.INd, n^+i — I.INd (a), no — I.INd, n^+i = 
3N D (b), n = 3N D , n N+1 = 1.1N D (c), and n = 1.2N D , n N+1 = 0.8N D (d) 



deviations for i = 0, 1 due to the contact are damped out because the field is in 
the PDC region there. In between there is a domain boundary where the electric 
field changes due to a charge accumulation e(n, — N&) = e r eo{Fi — -Fj_i). As 
the main jump occurs between Fig and F\q the density rij takes its maximum 
at ijj = 19. If now this charge accumulation layer is shifted by one period to 
%d = 20, we have almost the identical situation with the same current except for 
the fact that the voltage is diminished by (Fjjj — Fj)d as a lower fraction of the 
sample is in the high field region. This reveals the periodic sequence of branches. 
If one counts the branches starting from the right side as depicted in Fig. 0(a), 
the maximum of ni occurs in well id for the i£>th branch. In total one can 
count 40 branches which is exactly the number of quantum wells as the domain 
boundary may be located in each well. If the domain boundary comes close to 
i = l the low field region is not large enough to shield the variation due to the 
contact. Therefore the first and the second branch are strongly dependent on 
the boundary condition which simulates the contact (compare Fig. 0(a) with ||) . 
The field distribution in such stationary domain structures is not arbitrary as 
the currents across each barrier have to be equal. This provides a condition on 
the minimal doping density (or minimal carrier generation due to irradiation) 
as discussed in Refs. |54j, |6[ [|| |67| . 

Now the question arises why these stationary domain states are stable while 
in other spatially extended NDC systems like the Gunn diode typically travelling 
field domains occur. In Ref. 62] this question has been investigated using 
the simplified current relation (33). There it could be strictly proven that an 
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inhomogeneous field distribution is necessarily stable if all electric fields Fi are 
in the PDC region of v(F) which coincides with most readers' physical intuition. 
Therefore a field distribution like that in the inset of Fig. ^(a) must be stable if 
the NDC region eF max < eF < eF mnl is crossed within one jump. Then 

e-Fic-i < eF max and eF lD > eF min (34) 

holds. In the stationary state the current across each barrier has to be equal 
to Iq. Especially, there is Iq — I(Fi D ,rii D ,ni D +i). As Fi D is in the high field 
region the approximation ( |33| ) may be justified yielding I{Fi D1 rii D ,ni D+ i) w 
I(F iD , ND,ND)n iD /N£>. Combining this with Poisson's equation e{rn D —Nd) = 
e r eo(Fi D — F. tD ^i) exhibits that the condition (Q) can be fulfilled if 

Jo > I(F nlin ,N D , N D ) (l + £r£ ° (F ^~ Fmax) ) (35) 

holds. It has to be stated that this is only a sufficient condition for stability as 
there are stable domain states where one field is located within the NDC region. 
Now I can not be larger than J(J 1 max , Nr>, Njj) as the low field domain can not 
carry a larger current. Therefore there is a minimum doping density 

^ ereo(Jmin Jmax) Jnin 

^ Jmax Jmin 

above which stable domain states with an accumulation layer exist. From the 
I(Fi,N D ,N D ) relation from Fig. | Njgg » 1.2 x 10 10 /cm 2 is estimated which 
is much smaller than Nd = 1.5 x 10 11 /cm 2 . Therefore stable domain states are 
expected in accordance with the experimental and theoretical findings. Note 
that the condition (||^) depends strongly on J m ; n which itself is strongly affected 
by nonresonant transitions as shown in Fig. |[ 

Up to now domain structures have been discussed where the high field do- 
main is located at the receiving contact. But of course there is the other pos- 
sibility that the high field domain is at the injecting contact, i.e., Fi sa Fjh 
for i < id and F{ s» Fj for i > in- Then Poisson's equation yields a deple- 
tion region rii < No at the domain boundary Such domains have 
been observed experimentally for highly doped samples |3S|. Theoretically, 
such domains could be both obtained from a simplified model J63| and from 
the microscopic model |5q] . This is shown in Fig. |^ using the parameters of the 
sample from Ref. For the boundary conditions n (F 0l Nr>, n\) = 0.95Nd, 
nN+i(FN,ND,njsr) = 0.95-/Vd domain states are found, where the high field 
domain is located at the injecting contact as shown in the inset, while for 
n (F , N D ,ni) = 1.05N D , n N+ i(F N , N D ,n N ) = l.05N D the high field domain 
is located at the receiving contact like in Fig. [7| The domain branches them- 
selves look very similar. 

Regarding the stability of such domain structures the argument given above 
yields a minimum doping 

T\rdcp e r^o(Jmin Jmax) Jnax /'Q7'i 

crit ~ ~ 7 _ J . \ 6 ') 

c- i max 1 min 
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Figure 9: Calculated current-voltage characteristic of the sample from Ref. ]39j 
for voltage sweep-up for the boundary conditions no = n^+i = 0.95-/Vd (full 
line), no = njy+i = 1.05Nd (dotted line) 



above which stable domain states with a depletion layer can exist. For the 
sample of Ref. |Q the estimation gives N c ^ w 4.4 x 10 1:L /cm 2 which is three 
times larger than Njy. Thus, such domain states are not expected to be stable 
and therefore should not be observed. For comparison the sample used in |39| 
has N D = 8.75 x 10 n /cm 2 and from Fig. 3 of Ref. |§ N^g w 2.5 x 10 n /cm 2 is 
estimated in good agreement with the observation of stable domain structures 
with depletion layers. 

Note that the proof of stability essentially relies on the discreteness of the 
system. In a continuous model the NDC-region can not be crossed without 
having any fields within this region at least for a small spatial interval. This 
explains the difference to continuous systems like the Gunn diode where such 
stable domain states with an arbitrary position of the boundary are not ob- 
served. 



4 Imperfect Superlattices 

All the theoretical current-voltage characteristics shown up to now exhibited an 
almost regular series of branches whereas in typical experiments the lengths of 
the branches vary. It is straightforward to assume that this is caused by irreg- 
ularities in the real superlattice as nothing is perfect. But then the question 
about the nature of these irregularities arises. At first there are two different 



possibilities: Regarding a wafer as sketched in Fig. 10 the irregularity may ei- 
ther be a bad spot localized somewhere in the superlattice as shown in the right 
side of the wafer or a deviation from periodicity occurring in the whole layer. 
Such deviations may be a larger or smaller barrier width, or a different doping 
density in certain wells, e.g., which are established due to unsufficicnt control of 
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sample 1 sample 2 sample 3 

Figure 10: Schematic draw of a superlattice wafer from which 3 different samples 
are obtained. The dot in sample 3 represents a local imperfection while the line 
represents a modulation of the periodic structure. 



the growth process. Fortunately, one can distinguish these two cases experimen- 
tally. If the local spot would be the essential cause for the irregularities in the 
characteristics, sample 2 and sample 3 from the same wafer sketched in Fig. [To] 
should exhibit a different modulation of the branches as their individual spots 
have different sizes and are located at different positions. But the different sam- 
ples used in [Q exhibited almost identical modulations of the branches if they 
originated from the same wafers while there are large differences for samples 
from different wafers even if the superlattice structure is nominally identical. 
This would be expected from the samples 1 and 2 sketched in Fig. [h]. The 
same feature can be observed in |59) where the authors have fabricated samples 
with different numbers of periods originating from one wafer by an etching pro- 
cess. They found very similar sequences of longer and shorter branches which 
allowed them to conclude that the high field domain is located at the receiving 
contact in their samples. 

These experimental observations clearly indicate that the dominant devia- 
tions from periodicity are not (x, y)-dependent but extend over the whole wafer. 
This allows to simulate these irregular superlattices by introducing local fluctu- 
ations in the well width, barrier width or doping concentration into Eqs. (|l^ , p^ ). 
This has been done in Ref. [Q within the model of Ref. |33| where we found 
that even small spatial fluctuations (about 7%) of the doping have a significant 
influence on the length of the branches and can explain the observed behaviour. 
Furthermore, single fluctuations may be located by just determining the number 
of the branch which is altered most. This effect is particularly pronounced for 
fluctuations in the barrier width, where a one monolayer fluctuation may change 
the characteristic significantly. The theoretical prediction has been successfully 
checked experimentally by growing a new sample exhibiting one barrier with a 
larger thickness |68) . 

Now doping fluctuations are included into the model discussed before by 
using a well-dependent doping density Njji in Eq. (p9|). The impact on the 
impurity scattering is neglected, so that the old function I(F i; n.;, n i+ i) is used. 
The result is shown in Fig. 11 and exhibits fluctuating branch heights in ac- 
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Figure 11: Current- voltage characteristic for doping fluctuations, In (a) a ran- 
dom sequence Nm with an average fluctuation of 10% around No is used. In 
(b) there is -/Vdi = l.2Nn and Nd3o = 0.8Nd while the other densities are not 
altered. 



cordance with the model used in Ref. |Q. Comparing the local fluctuations 
Noi with the maximum current of the branch, one finds a correlation. Like the 
findings of Ref. |39) the branch i + 1 counted from the right extends to higher 
currents if Nm is larger than Hp. Additionally the branch i extends to lower 
currents as can be seen from Fig. |ll|(b). The deviations in the current seem 
to be smaller and the effects to neighboured branches seem to be larger here 
than observed in Refs. |6j, |69| . This might be related to the fact that a higher 
doping density was used in previous works. 
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5 Oscillatory behaviour 

Oscillatory behaviour has been observed in coupled multiple quantum wells 
both experimentally J26[ p7| , |28{ and theoretically 30 1. As an example such 
behaviour is obtained within the model discussed here for the boundary con- 
ditions no(F,Nn,ni) — 0.8Nd and tin + i(F 7 Nd, tin) — 1.2Nn which is just 
the reversed sequence used for the calculation of Fig. ||(d) where stable domain 
branches were found. Fig. [l2| shows self-sustained current oscillations between 
7.4 /iA and 8.2 /iA with a frequency of 0.18 MHz. Similar results are obtained 
for different biases and different values no(F,Nr>,ni) < No- The oscillations 
may be described from the field distribution as follows. At t = 5.71/xs the field 
distribution consists of a high field and a low field domain with a depletion layer 
in between. As for Ne> < N c ^ such a distribution is not stable the boundary 
travels to the right thereby decreasing the electric field in the high field domain 
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because the total voltage has to remain constant (t — 8.45^s). As the high 
field domain enters now the region of NDC for eFd < 90 meV a new positive 
charge accumulation layer is created there which is slightly visible for t = 9.21/is 
at i « 10. This accumulation layer travels to the right and increases in time 
(t = 10.16/is) until it merges with the old depletion layer. Then the cycle is 
repeated again. Such an oscillation type has been described in Ref. |n) within 
the approximation (p3|). The very same behaviour happens for charge accumu- 
lation layers as well if domain states with an accumulation layer are unstable. 
For a full treatment of this oscillation type see Ref. |32|. Note that the mini- 
mum current in this oscillation cycle roughly coincides with I(F mm , Njj, Njj). 
This relation seems to hold generally as checked by altering the nonresonant 
current and thereby I(F m i n , Njj, Njj). Therefore the temporal minimum of this 
type of current oscillation provides information about I(F mm , No, Nd) which 
is strongly dependent on the nonresonant transitions. 

In order to study the influence of the boundary conditions onto the oscilla- 
tions we use the ohmic contact currents /o^i = ueF^d and In^n+i = ueF^d 
in the following. For a = 0.5 mA/eV a completely different oscillation mode 
is found as shown in Fig. |l3|. Here one oscillation cycle consists of the nucle- 
ation and travelling of a dipole domain which vanishes by leaving the sample at 
the receiving contact. The scenario is completely analogous to that described 
in |?(J for a continuous system exhibiting NNDC. Additionally, small current 
spikes appear which are related to the motion of the accumulation layer from 
one well to the next as discussed in Ref. ]7lJ for the domain formation process. 
Thus, these spikes reflect the discreteness of the superlattice. 

Experimentally, the sample discussed here also exhibits self-sustained current 
oscillations if a positive bias is applied at the top contact while there are stable 
domain branches for a negative bias. The experimental data j3lj exhibit current 
oscillations between 20 /j,A and 35 fiA with a frequency of 0.4 MHz at U = 0.78 
V. Thus, the theory is in qualitative agreement with the experiment regarding 
the oscillations as well. As we did not try and model the contact currents 
microscopically (which strongly influence the type of oscillations) we can not 
expect quantitative agreement. Additionally, the shape of the oscillations as 
well as the frequency depends substantially on the full shape of the current-field 
relation which is strongly affected by the nonresonant transitions between the 
resonances. Furthermore, deviations from periodicity may affect the oscillations 
as well @ [[§. 

For larger values of a (say a > 3 mA/eV) the same domain branches like 
those shown in Fig. ^ are found. This indicates that the domain branches them- 
selves are almost identical if their formation is allowed for by the contact con- 
ditions and if the doping is sufficiently large. On the other hand the oscillatory 
behaviour does strongly depend on the exact contact conditions. 
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Figure 13: Calculated current oscillations for U = 1 V using ohmic boundary 
conditions with a = 0.5 mA/eV. 

6 Details of the calculations 

In this last section I want to present the details of the calculations whose results 
have been shown before. 



6.1 Calculation of the couplings 

In order to calculate the coefficients of Table |l| for Eqs. (13-15]) we have to 
specify the band structure in GaAs and AlAs at first. It is assumed that only 
the r band is of importance. While one may use a parabolic band structure 
for ALcGai^As/GaAs heterostructures for small x this is not appropriate for 
an GaAs/ AlAs heterostructure as the conduction band of GaAs is located far 
in the band gap of AlAs where the band structure is clearly not parabolic 
(see Ref. [fnj). Following Ref. [M we model the nonparabolicity by an energy 
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dependent effective mass m(E) = m c (l + (E—E c )/E g ), where m c is the effective 
mass at the conduction band minimum of energy E c , and E g is the energy gap. 
Then the usual connection rules 

tp(zo-e) = tp(z + e) (38) 
1 d ^o-e) = * , A zo + e) 09) 



m(E, z — e) dz m(E, z + e) dz 

hold for the envelope function (p(z) provided that the momentum matrix element 
P = h-\f E g /(2m c ) between the conduction and valence band states is identical 
in both materials. We use the values 173] m[f aAs = 0.067m e , m A1As = 0.15m e , 
^GaAs _ 52 eV, iJ A1As = 3.13 eV, and the conduction band discontinuity 
AE C — 1.06 eV. These parameters yield a relation E(k) = E c + h 2 k 2 /(2m(E)) 
which is in excellent agreement with the band structure of AlAsj73| for the en- 
ergies of interest. It must be stated that for these parameters the value of P 
is slightly different for the two materials in contrast to the assumption. Never- 
theless, these parameters seem to give reasonable agreement with experimental 
data both in Ref. Jt3] and the calculations presented here. 

Imposing the Bloch condition (fg(z + d) = e lqd (p g (z) the Bloch functions 
<Pq(z) and eigenvalues E v q are calculated within the Kronig-Penney model. The 
phase of the Bloch functions is chosen in the following way Q : Let z — be the 
center of one quantum well. If ^ we chose the phase in such a way that 

</j^(0) is real for each q. For <^q(0) = 0, <Pg(0) is chosen to be purely imaginary. 
Furthermore ^p v q (z) has to be an analytic function in q for both cases. From 
Eq. the level energies E v and couplings T" are obtained. Eq. (|l|) provides 
the Wannier functions which are plotted in Fig. ^. Finally, the couplings i?^ " are 
obtained from their definitions W h v — J dz^ v (z — hd)z 1 $ l/ (z). The calculated 
values are given in Table [l]. A complication arises due to the fact that the 
effective Hamiltonian of the Kronig-Penney model is energy dependent due to 
the energy dependence of the effective mass. Therefore the envelope functions 
(fq(z) for different energies E^ are not strictly orthogonal but exhibit a small 
overlap which is neglected in the calculation. 

6.2 Impurity scattering 

Here we want to calculate the self-energy 5 for impurity scattering. The contri- 
bution to H* ca - tteT of Eq. (|l^) for the lowest level a is given by 

^ P =J E V^(p)ai(k+p)a n (k) (40) 

Here the subscript i denotes the impurity located at the position (r\j, z{). The 
matrix element is calculated with the Wannier functions yielding: 

„2 



V%(p) = J d 2 rdze~^- K{z) y a{z) 



4TT€ s e ^J\r-r % \ 2 + (z - Zi + nd) 2 



5 For readers who are not familiar with the concepts of many particle physics (such as 
self-energies, Green functions, etc.) Ref. |75| is recommended as a helpful introduction. 
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n\ V V <.n+h n* in 



/ TT tULL' V'V 
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(b) 



v it 

rffi v r' 






n' iK—^ai' riv_^n n+hV Vn' 



Figure 14: a) gives the notation of the matrix element W£ (p) and the 

polarizability 11%"$ " {p,u = 0). b) shows the random-phase approximation 
diagramatically 

_ p 2 /• 

= —— / dz^* a {z)^ a (z)e-^ z - z ' +nd ^e- l P^. (41) 

Note that the contribution for p = is canceled by the respective part in the 
electron-electron interaction as usual if the number of donors is equal to the 
number of carriers in the whole sample. In order to prevent the divergence of 
the matrix element for p — > screening due to the electron-electron interaction 
has to be considered. (This is a general problem for the Coulomb interaction, 
see also chapters 5 and 6.) The Hamiltonian for the electron-electron interaction 
reads: 

H ee = -jj Wr a i^<(k + P)< +h {}l~PW + H(kL)a n {k) 

k.k' ,y,n,h 

+terms with b . (42) 

Here we assumed that the overlap of wave functions from different wells is 
negligible, so that pairs always have to be inside the same well. The matrix 
element reads: 

WT V » = 7T~~~ — / dzr [ dz 2 ^*( Zl )^'*(z 2 - dh) 



2e s e Q p 



x^'(z 2 ~dh)^(z 1 )e- plzi ~ Z21 (43) 



The screening is described by the polarizability v (p, ui = 0) , where v, \x 

take the values a and b. The notation follows Ref. pq , where a similar problem 
is investigated, and is shown in Fig. [l4|. As the p dependence is identical in 
all parts, we omit it in the notation. Within the random-phase approximation 
(RPA) we have the Dyson equation (see, e.g., Rcf. ]5q| ): 

jjHV, fj,' u' jjO/iv I r s r , yrun" u" Hjrii" u" ,n' v' \ 

ll n,n' — Ll n \ °n,n' li,,v'Ou,n' T 2_^i VV h 11 n+h,ri I ' 

h,p",v" 
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where the vacuum polarizability II is given by pSj 

n - A k {E k+R + E») - (E k + E») V 7 

Now we find that IT° bb (p,w = 0) = assuming that the upper level is not 
occupied. II° afc (p, u = 0) is quite small as the gap appears in the denominator. 
It remains the contribution IT° aa (p,a; = 0) which yields for the 2 dimensional 
electron gas Q : 



n°"(p, W = 0)=-^ I 1 - Q(p - 2fcp)y 1 - 4^ | (46) 

which is independent of p for p < 2kp. As all polarizations only have a- indices 
we omit these indices in the following. Eq. (0) can be solved by a Fourier trans- 
formation (for an infinite superlattice and assuming 11° = II is independent of 
n). Defining IL q = J2 n > H n , n 'e iq( - n '~ n) , W^ aaa = J2h W£ aaa e lqh we find 

n° 

n, = = . (47) 

Now the screened electron-impurity interaction is given by the bare interaction 
and a part combined with the polarizability given by 

V™ sc (P) = Vg + Wr a V- n+h ,n' . (48) 

h,n' 

Defining V?® = ^2 n V^ n e~ l9 ( n ~ ni ) and using the translational invariance H n +h,n' 
Tl2n+h-n',n we find: 



yaa 

i,<l iq(n—rii) 



(49) 



Similarly we have: 

v? n sc ( E ) = v% + ]T wj; aab n n+h , n ,v™, (50) 

n' 

~\fbb /ytrbaabTyaa ^raaaa~Crbb\T~\^ 



rbb iTxrhaab^raa ytraaaa^rbb\y\0 

' i,q . \ q i,q q i,Q 



iqin—rii 



The screened matrix elements from Eqs. ( (f9|j5C| ) are used in the Hamiltonian ( [40[ ) 
in the following. 

Now we calculate the self energy within the single-site approximation which 
is shown diagrammatically in Fig. nH. For the impurity i we find the contribution 



2N 



£ = / \+ / l\ + // \\ +... 

' \ / I \ , ' \ v 
( ^_ i t ^ i ^ \ ' ^ L_^_ i ^ N 

— ^ — — 7^ ^ — — ^ ^ 7^ — 

Figure 15: Diagrams contained in the self-consistent single-site approximation 
to the self energy S» ret (fc, E) 

K'Hk, E) = j^Y, V *T C & ~ hMk^E^rik, - k) (5i) 

+js E ^r c (fc-fci)G(fc 1 ,s)^ sc (fci~fc 2 )G(fc 2 ,i?)^ sc (fc2-fc) 

+ Ji E y i T C (k-k 1 )G{k 1 ,E)V™ sc {k 1 ~k 2 )G{k 2 ,E) 

x^r c fe - fc 3 )G(fc 3 , s)^ sc (fc 3 - *) + . . . 

where G(k, E) = (E - E k - S^ rct (A, is the full retarded Green function. 

This sum can be transformed to the self-consistent equation (see, e.g., [(59|) 

K^ih, k, E) = V™ sc (ki ~ K) + j E Kn "Ui - k 2 )G(k 2 ,E)K a '%k 2 ,k, E) 

(52) 

which can be solved numerically for a given self-energy function S° ret (fc, E 1 )) 
entering G{k 2 , E). The contribution to the self energy is then given by 

K' l (k,E) = ^J2 V ^ SC (h-k)G(k 1 ,E)K a '% 1 ,k,E). (53) 

Summing up the contribution from all impurities i and possibly different scat- 
tering processes we obtain the self-energy function 

K ICt (k, E)=J2 K'Hk, E) + Smother scattering ^ ^ 



Now a self consistent solution of Eqs. (52, 53] , |54| ) can be achieved by iteration. 
The contribution to S^ ret (fc, E)) is calculated in the same way. 

6.3 Interface roughness 

Interface roughness is modelled like in Ref. J56| considering an interface located 
at z = zo exhibiting thickness fluctuations £(r) of the order of ±77 (we use 
77 = 2.8 A which is one monolayer of GaAs). We assume the correlations 

(e(r)>r = (55) 
<£(r)£(r')>r = a^expHr-r'IM) (56) 
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with a correlation length A = 7 nm and an average coverage a = 0.5. Such an 
exponential correlation function (£(n)£(n'))r seems to be more appropriate than 
the usual choice of a Gaussian (see, e.g., Ref. |78]|) as stated in Refs. |79, 56). Like 



in Ref. [J78l (where the scattering of Bloch states in a superlattice is regarded) 
the additional potential^] due to the roughness is modelled by a ^-function at 
the perfect interface 

U(r, z) = i(r)AE c 6(z - z ) . (57) 
This gives the following interface roughness contribution to jy scatter 

trough = ^ [^(^o|, +h (*+£)o„(^+Cr»(£)6t +h (fe+g)6 n (>0 
k,p,h 

+U b h a {p)bi +h {k + p)a n (k) + Ufi b (p)ai +h (k + p)b n (k)] (58) 
with the matrix elements 

U7(P) = / d 2 re' l P- r -AE c [e(r)*"*(z - M)*^o)] • (59) 



Using the correlation function ( |56"| ) we obtain the square of the matrix element 

\U^(p)\ 2 = AAE*\*„(zo - hd)| 2 |*M*o)| 2 , W ^ 3/2 (60) 

(1 + (pA) ) 

which enters the expressions in the following. 

The elements £/q M result in scattering within the wells. Their contribution 
to the self-energy is calculated within the self-consistent Born approximation 
(which is just the first diagram of the infinite sum in Fig. [l5| ) 

£ r a r u c gh(fc, E) = \ U ° a & - t^Gik^E) , (61) 

where the factor 2 takes into account the two interfaces per well. These self- 
energies contribute to the total self energy in Eq. (|5~4|). The calculation for the 
subband b is performed in the same way. 

The elements IL contribute to the non-resonant current from one well to 
the next via Eq. ( J18| ) . Here the contributions from all 4 interfaces of both wells 
involved are summed up. For weakly coupled wells , J7g M , . . . are small and 
can be neglected. 



6 In chapter 11 of this book the potential is chosen to be the variation of the energy levels 
due to the well width fluctuation which is only defined for scattering within a given level 
and a given well. The approach via Eq. has the advantage that interwell and interlevel 
transitions can be handled as well. For intrawcll and intralevel processes the results are similar 
as AE V /dw ~ A_E C | \P" (zq ) 1 2 where w is the well width and zq is the position of the interface. 
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6.4 Optical phonons 

In polar materials like GaAs the polar interaction with optical phonons provides 
an important scattering process. As the energy huj — 36 meV is transferred, 
this process couples electronic states with different energy in contrast to the two 
scattering mechanism discussed before. This makes the full calculation of the 
self-energies much more complicated. But fortunately only a restricted number 
of processes are allowed at low temperatures where only the emission of phonons 
takes place. Therefore to any electronic state k in level v with energy Ek + E v 
affected there must be another state at an energy Ek + E v — Hujq. For the lowest 
level v = a this means that the condition Ek > hcuo must be satisfied in order 
for phonon scattering to be possible. But for small Fermi levels (5.37 meV for 
the sample considered) these states are far away from any resonant transition so 
that the neglection of the phonon contribution to the self-energy hardly affects 
the currents. The situation is different for the second level. Here the states 
with Ek ~ are in resonance with the occupied states in the ground level if the 
electric field takes the value eFd as E b — E a . Therefore the actual broadening 
of these states is crucial for the a — > b resonance. This process is taken into 
account by calculating the scattering time r p h for this process following Ref. |pl[ 
yielding r p h = 0.854 psec for the structure considered here. For the self-energy 
contribution the constant value 

SjhL(fe^)=-i^- (62) 
is used which contributes to the total self energy in Eq. (E3) for the level b. 

7 Conclusions 

In this chapter the electronic transport in weakly coupled multiple quantum 
wells has been considered. Within the model of sequential tunnelling the cur- 
rents have been calculated without any fitting parameters taking into account 
the scattering at impurities and interface roughness. The currents are in good 
quantitative agreement with experimental data stating the physical relevance of 
the model. 

In the NDC region a homogeneous field distribution is unstable. For the ac- 
tual doping of the sample considered both stable field-domains and self-sustained 
current oscillations are found theoretically in good agreement with the exper- 
imental observation. The sequence of domain branches is almost independent 
of the contacts as the influence of the boundaries is shielded by the domain 
regions where the electric field is in the regime of positive differential conduc- 
tivity. Therefore the branches contain information about the transport in mul- 
tiple quantum wells itself which is not spoiled by contacts which are often only 
poorly defined. Furthermore the domain branches react quite sensitively to lo- 
cal deviations from periodicity which allows for a check of the actual sample 
quality. 
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